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Abstract 

In the previous paper we studied the transport coefficients of Quark-Gluon- 
Plasma in finite temperature and finite density in vector and tensor modes. In this 
paper, we extend it to the scalar modes. We work out the decoupling problem and 
hydrodynamic analysis for the sound mode in charged AdS black hole and calcu- 
late the sound velocity the charge susceptibility and the electrical conductivity. We 
find that Einstein relation among the conductivity the diffusion constant and the 
susceptibility holds exactly 



1 Introduction 



The discovery of low viscosity in the theory with gravity dual [1] and its possible relation 
to the RHIC (Relativistic Heavy Ion Collider) experiment induced a great deal of efforts to 
establish the relevant calculational scheme that may be provided by AdS / CFT correspon- 
dence [2-4]. An attempt has been made to map the entire process of RHIC experiment in 
terms of the gravity dual [5]. The way to include a chemical potential in the theory was 
figured out in the context of probe brane embedding [6-13]. Phases of these theories were 
discussed and new phases were reported where instability due to the strong attraction is a 
feature [8-10]. 

In spite of the difference between QCD and M = 4 SYM, it is expected that some 
of the properties are shared by the two theories based on the universality of low energy 
physics. In this respect, the hydrodynamic limit is particularly interesting since such limit 
can be shared by many theories. The calculation scheme for transport coefficients is to 
use Kubo formula, which gives a relation to the low energy limit of Wightman Green 
functions. In AdS/CFT correspondence, one calculate the retarded Green function which 
is related to the Wightman function by fluctuation-dissipation theorem. Such scheme has 
been developed in a series of papers [14-18]. 

For the hydrodynamic analysis, one may need to have master equations for the decou- 
pled modes in vector and scalar at hands. Although the analysis for the decoupling problem 
were analyzed in [19], it was based on 5*0(3) decomposition while more useful work for 
hydrodynamics should be based on 5*0(2) decomposition, where longitudinal direction of 
the spatial direction is distinguished. For this purpose, Kovtun and Starinets worked out 
the decoupling problem based on 50(2) for the AdS black hole case [18] before doing the 
hydrodynamic analysis. For the charged cases, there are extra difficulties: vector modes 
of gravity and those of the gauge fields couple. Furthermore there are extra couplings in 
scalar modes which are not present in the chargeless cases. 

In the previous paper [20], some of us extended this work to charged case using the 
Reissner-Nordstrom-Anti-deSitter (RN-AdS) black hole, which corresponds to the diagonal 
(1, 1, 1) R-charged STU black holej. However, analysis for the scalar mode of charged case 



* In fact much works had been done for charged case by various groups [21-25]. In [21,22], ther- 
modynamics for STU black hole [26,27] and the hydrodynamic calculations for the (1,0,0) charge were 
performed. In [23,24], charged AdSs and AdS4 black hole backgrounds were considered, respectively, and 
it was shown numerically that the ratio (rj/s) was l/(4.7r) with very good accuracy. Later, it was also 
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was not done due to difficulties caused by extra mixing between various scalar modes in 
charged AdS black hole. In this paper, we work out the decoupling problem and hydrody- 
namics for the sound (scalar) mode in such case. Green functions are explicitly obtained. 
Our results show that the behavior of the transport coefficients in RN-black hole are very 
different from those in the (1,0,0) charged black hole: the formers are much more smoother 
than the latters. We find that Einstein relation among the conductivity, the diffusion 
constant and the susceptibility holds exactly. 

This paper is organized as follows: In Section 2, we introduce RN-AdS black hole 
and review correlation function calculation at finite temperature in AdS/CFT correspon- 
dence. In Section 3, a formulation on the metric and the gauge perturbations in RN-AdS 
background is considered. We then solve linearized perturbative equations of motion in hy- 
drodynamic regime and obtain retarded Green functions in Section 4. We also observe the 
transport coefficients including the speed of sound, the diffusion constant for U(l) charge 
and the electrical conductivity in this section. Conclusions and discussions are given in the 
final section. Three appendices are provided. In Appendix A, we summarize the results 
in the vector and the tensor type perturbations in our previous work [20]. The details of 
calculations to solve equations of motion are given in Appendix B and C. 

2 Basic Setup 

2.1 Minkowskian correlators in AdS/CFT correspondence 

Before introducing RN-AdS black hole, we briefly summarize to obtain Minkowskian cor- 
relators in AdS/CFT correspondence. We follow the prescription proposed in [14]. We 
work on the five- dimensional background, 



where x M and u are the four-dimensional and the radial coordinates, respectively. We refer 
the boundary as u = and the horizon as u — 1. A solution of the equation of motion 
may be given, 



ds 2 = g^dx^dx" + g uu (du) 2 , 



(2.1) 




(2.2) 



proven that the ratio might be universal in more general setup [25]. 
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where fk(u) is normalized such that /^(O) = 1 at the boundary. An on-shell action might 
be reduced to surface terms by using the equation of motion, 

J {^) u=0 

Here, the function Q(k,u) can be written in terms of f±k(u) and d u f±k{u). Accommo- 
dating Gubser-Klebanov-Polyakov/Witten relation [3,4] to Minkowski spacetime, Son and 
Starinets proposed the formula to get retarded Green functions, 

G R (k) = 2g(k,u) , (2.4) 

where the incoming boundary condition at the horizon is imposed. In general, there are 
several fields in the model. We write the Green function as Gij(k), where indices % and j 
distinguish corresponding fields. 

In this paper, we work in RN-AdS background and consider its perturbations so that 
essential ingredients are perturbed metric field and U(l) gauge field. Here we define the 
precise form of the retarded Green functions which we discuss later: 

pa {k) = -i J & A x e- ikx 9(t){[T^(x), T PCT (0)]>, 
G, v p (k) = -ij d 4 x e- ik *9(t)([T^(x), J„(0)]>, (2.5) 
G, v {k) = -i J & A x e-^e(t)([J,(x), J„(0)]), 

where the operators T^ u (x) and J M (x) are energy-momentum tensor and U(l) current which 
couple to the metric and the gauge field, respectively. 

2.2 Reissner-Nordstrom-AdS background 

The charge in RN-AdS black hole is usually regarded as i?-charge of SUSY [28]. We here 
consider an another interpretation in the following way: One can introduce quarks and 
mesons by considering the bulk-filling branes in AdS 5 space. The overall U(l) of the flavor 
branes is identified as the baryon charge. The U(l) charge in this model [29] minimally 
couples to the bulk gravity since the bulk and the world volume of brane are identified. 
Then, the baryon charge and the i?-charge have the same description in terms of the U(l) 
gauge field living in the AdS 5 space. A charged black hole (RN-AdS black hole) is then 
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induced by its back reaction. Therefore the U(l) charge in RN-AdS can be identical to the 
baryon charge. As a result, we can interpret our result as a calculation of the transport 
coefficients in the presence of the baryon density. 

The effective action of this gauge field is given the quadratic piece of Dirac-Born-Infeld 
actior * 



^ J d 5 x^ Tr {T mn T mn ), (2.6) 
where the gauge coupling constant e is given by [29] 



Sdbi ~ Ae 2 



I N c N f , 

? = -U< (2 - 7) 

with I the radius of the AdS space. We pick up an overall U(l) part of this gauge field in 
order to consider the baryon current at the boundary. Together with the gravitation part, 
we arrive at the following action which is our starting point: 

S ° = ^l d " X ^ { R ~ 2A ) " hi ^V^Q^mn^, (2.8) 

where we denote the gravitation constant and the cosmological constant as k 2 = 8ttG^ and 
A, respectively. The U(l) gauge field strength is given by J r m n{ x ) — d m A n (x) — d n A m (x). 
The gravitation constant is related to the gauge theory quantities by 

Suppose we have baryon charge Q. This should be identified to the source of U(l) charge 
on the brane hence on the bulk. Then we can relate it to the parameter in RN black hole 
solution by considering the full solution to the equation of motion, 

Rmn ~^QmnR ~\~ Qmn^- ^ T mn ^ (2.10) 

where energy-momentum tensor T mn (x) is given by 

An equation of motion for the gauge field A m (x) gives Maxwell equation, 

W m r mn = -^d m (y^g mh g n \d k Ai - diA k )) = 0. (2.12) 



* The indices m and n run through five-dimensional spacetime while /x and v would be reserved for 
four-dimensional Minkowski spacetime. Their spatial coordinates are labeled by i and j. 



Here we assumed that there is no electromagnetic source outside the black hole. One can 
confirm that the following metric and gauge potential satisfy the equations of motion (I2.10p 
and ([ZZEZD , 

d* 2 = £( - f(r)(dtf + E( d ^) 2 ) + ^)( dr ) 2 ' ( 2 - 13a ) 



i=i 



A t = ~+fi, (2.13b) 



with 



if and only if g is related to the Q by 

« 2 - |p- (2-14) 

It should be noted that a ratio of the gauge coupling constant e 2 to the gravitation constant 
k 2 is 

e 2 iV c 9 , 

- = *r' ■ 2 - 15 

«r TV/ 

Since the gauge potential At(x) must vanish at the horizon, the charge Q and the chemical 
potential // are related^. The parameters m and q are the mass and charge of AdS space, 
respectively. This is nothing but Reissner-Nordstrom-Anti-deSitter (RN-AdS) background 
in which we are interested throughout this paper. 

The horizons of RN-AdS black hole are located at the zero for /(rjl, 

ml2 , 1 ( 2 2\ ( 2 2 

- -r + )[r -r_ 

where their explicit forms of the horizon radiuses are given by 



f(r) = 1 - ^_ + 1_ = _ ^ _ r% j ^ _ r ^ ) ( r * _ r g ) , (2 . 16) 



r + = (wA 1+2cos ^ + r))) • (2 - 17a) 

r 2 = (^(l + 2cos(t) ) ) , (2.17b) 



3g2 ^- ■ V3 

773 

V3g 



-o=fefl + 2cos^ + ^))) . (2.17c) 



' The chemical potential /i can be expressed by using gauge invariant quantity as 

poo 

fJ,= ArT r t = At(aa), 

J r+ 

where r + and oo represent the horizon and the boundary, respectively. This definition gives thermodynamic 
relations consistently. 

* In order to define the horizon, the charge q must satisfy a relation q 4 < 4m 3 l 2 /27. 



with 

V 2m 3 / 2 - 27g 4 ) ' 
and satisfy a relation r\ + r 2 = — r 2 . The positions expressed by r + and r_ correspond to 
the outer and the inner horizon, respectively. It is useful to notice that the charge q can 
be expressed in terms of 9 and m by 

4 4m 3 / 2 . 2 



9 = ^ sin - 
H 27 V 2 

The outer horizon takes a value in 

/ jyi 

J—l<r 2 + < y/ml, 

where the upper bound and the lower bound correspond to the case for q = and the 
extremal case, respectively. 

We shall give various thermodynamic quantities of RN-AdS black hole [28,29]. The 
temperature is defined from the conical singularity free condition around the horizon r + , 



rlf'(r + ) _ rw 1 q 2 i 2 \ _ 1 



r = J sr i = ^ 1 -27rJ s 2sl 1 -2J- (>0) - (2 - 18) 



a 



where we defined the parameters a and b by 

b = 

2r+ 



a 2 l 2 I 2 
a=^r, b=^-. (2.19) 



' + 

In the limit q — > 0, these parameters go to 

a — > 0, 

and the temperature becomes 

T ^ T °-nP/ 2 - 

The entropy density s, the energy density e, the pressure p, the chemical potential p and 
the density of physical charge p can be also computed as 



I 3 / 2 
^ 2mV4' 

m l/4 





2vrr 3 


tt/ 3 


s = 


k 2 / 3 


46 3 /t 2 ' 




3m 


3/ 3 


e = 


2/t 2 / 3 


326 4 k 2 




e 




P = 


3' 






Q 


46 2 Q 


/i = 


r 2 " 


Z 4 ' 




2Q 


/ p 


P = 


e 2 / 3 


= ^2b^' 



(l + a), 



(2.20) 

(2.21) 
(2.22) 

(2.23) 
(2.24) 



In order to obtain a well-defined boundary term from the gravitational part, we have 
to add the Gibbons- Hawking term into the action, which is given by 

S ^ = ^j d A xy/^K, (2.25) 

where integration is taken on the boundary of the AdS space. The four-dimensional metric 
9$(x) is the induced metric on the boundary and K(x) is the extrinsic curvature. We also 
need to add counter terms to regularize the action [30], 

5 ct = ^/d 4 xv^)(y-^). (2.26) 

3 Perturbations in RN-AdS Background 

In RN-AdS background, we study small perturbations of the metric g mn {x) and the gauge 
field Am(x), 

a = o (0) + h 

iJmn — <Jmn "-rrm; 



A = A®) j- A 



(3.1) 
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where the background metric gmn{x) and the background gauge field Am\x) are given in 
f l2.13al) and (12.13bj) . respectively. In the metric perturbation, one can define an inverse 
metric as 

g mn = g (0)mn _ frmn + ^nl^n + Q(h 3 ), 

and raise and lower indices by using the background metric gmli(x) and g^ mn (x). 

Let us now consider a linearized theory of the symmetric tensor field h mn (x) and the 
vector field A m (x) propagating in RN-AdS background. We shall work in the h rm (x) = 
and A r (x) = gauges and use the Fourier decomposition 

d A k 

d 4 k 

where we choose the momenta which are along the ^-direction. In this case, one can 
categorize the metric perturbations to the following three types by using the spin under 
the SO (2) rotation in (x, y)-plane [15]: 

• vector type: h tx ^ 0, h zx ^ 0, (others) = 

^equivalently, h ty ^ 0, h zy ^ 0, (others) = 



MW) = / 7^1 e- Mz h, u (k,r) : 



A„(t,z,r)= e ~ Mkz A,(k,r) 



• tensor type: h xy ^ 0, (others) = 

^equivalent ly, h xx = —h yy ^ 0, (others) = 

• scalar type: h tt ^ 0, h tz ^ 0, h xx = h yy ^ 0, and h zz ^ 0, (others) = 

First two types of the perturbations were studied in [20]. We list the result in Appendix 
A. In this paper we consider the scalar type perturbation. 

3.1 Linearized equations of motion 

From explicit calculation, one can show that t and ^-components of the gauge field A^(x) 
could participate in the linearized perturbative equations of motion. Thus independent 
variables are 

h tt (x), h tz (x), h xx (x) = hyy(x), h zz (x), 
A t (x), A z (x). 

In the hydrodynamic regime, it is standard to introduce new dimensionless coordinate 
u = r+/r 2 which is normalized by the outer horizon. In this coordinate system, the horizon 
and the boundary are located at u — 1 and u — 0, respectively. We also define new field 
variables 

h t _ n (o)tt h _ l 2u h 
n t — 9 n tt — Tl n tti 

nf 

K = g (0)zz h zt = l ^h zt , 

h x = J0)xx h = L _^ h 

,u x y ,b xx o ,b xxi 

If = a (0)zz h = l ^_ h 
'"z y "-zz 2 

r + 

A I 4 
r — z_itL — a 

fx 

where \x is the chemical potential given by (12.231) . Nontrivial equations in the Einstein 
equation (12. 8p appear from (t,t), (t,u), (t,z), (u,u), (u,z), (x,x) and (z,z) components, 
respectively: 

u-ft t + u _ lf [ 2 n t +n x + 2 n z ) uf ti t 

f — [u 2 h x x + -u 2 h z z + ukh{) + 2a-h\ + AajB' t , (3.2a) 
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o = u Uhf + K'-j (K + \K) \ + k (hf - f jK 

= hf - -hf + =^hZ - SauBL 



u 



= hf + 2hf + hf + f j (hf + hf + \hf) + 2a U - f h\ + 4ajB' t 
= khi + 2khf - jhf + jjrhl + 3a j (kB t + uB 2 



= h x " + ^—ILh x ' ■ 

U ll X ' or ' l X 

U 2 f 

-a-h\-2a-B' t , 
f f 

= hf + ^lh? 

U 2/ 

<x—h\ - 2a-B' t , 
f f 



2u\ 



( hi' + hf 



+ 



uf' 



(3.2b) 
(3.2c) 
(3.2d) 
(3.2e) 

(3.2f) 



u\2 X J uf- 



(u 2 h z z + 2ukh z t - k 2 fh\ - 2k 2 fh x x 



(3-2g) 



with 



f(u) = (1 -u){l + u-au 2 ), 

where the prime implies the derivative with respect to u. On the other hand, in the Maxwell 
equation (12.121) . t, u and ^-components give nontrivial contributions 

= B'i - (k 2 B t + kuB z ) + i (hf - 2hf - hf) , 



= uB> + kfB' z + ^(hl- 2hl - hi) - khf 
= B 'i + f -B' z + (uj 2 B z + tjkBt) - jhf. 



(3.3a) 
(3.3b) 
(3.3c) 



The equations (I3.3al) and ( I3.3bl) imply (I3.3cp . In the set of equations for the metric per- 
turbation ( 13.2al) -( |3~2"g"| ), together with ( 13.3al) and ( I3.3b[) . the following four independent 
relations are obtained: 



3(u 2 - k 2 f) + k 2 uf 
uj( y f'(3f-uf')-Ab 2 LU 2 



h ' 2k\3f-uf>) 



2b 2 u 2 hX _ f'(3f-uf>)-4bW ht 



4/(3/ -uf) 



2k 2 P(3f-uf) 
3au 2 u 2 (. „, 



+ -K + kh z t 



3au 



2k*f(3f-uf) 

hZ , = 3lo 2 + h*uf ht , 2b 2 P 



k*(3f-uf) 1 3f-uf 



2kf 



kB t + uB z j , 
2b 2 



(3.4a) 



-(c 2 (2^ + ^)+2^ 



-^^y ((/W - uf) - M)* - W K ) 
cu(f(3f-uf)-4b 2 u; 2 ) 



k 2 f 2 (3f-uf>) 
3au 2 (u 2 + k 2 f) (ut a , 3a M 



_ 3^/ , 2i\j / N /' (3/ - «/') - 4iV / u 

'' ~ k(3f- uf) ' + 3/ - «/' I ' + K ) + kf(3 f - uf) \ * + 2 ' + kh < 

W -( ft S + 2B ;), (3.4c) 



k(3f - uf) 



O - ft '' + 2 M /(3/- M /0 ^ 3(/ "" /,)(2/ "" m ' 



+46 2 ^ - fc 2 //i* + (2u 2 + (f- uf')k 2 ^ h x x + u 2 h z z + 2u;kh^j 
+au 2 (15 f - 7uf) fa + }. (3.4d) 



The equations of motion (l3.2al)-(|3T2"g"D can be derived by using the above relations. Taking 



the limit q — > 0, the relations ( 13.4aP -( l3.4dl) coincide with the result in [18]. 
3.2 Surface terms 

Before solving the equations of motion, we shall give a surface action in oder to obtain 
Green functions. By using the equations of motion, bilinear parts of on-shell action (I2.8P 
reduce to surface terms: 

I 3 f d 4 k ( f t t , f , f , 3 , 



s ° = Wva / T^li +-W + -KK + -KK - -Mh 

32k 2 o 4 J (27r) 4 I u u u u 

-Lh x h 1 ' —tfh 1 ' - —h t h xl 

u 2u u 

-Lh z h x ' —tfh*' - —h x h zl 

,L z' L x o u t u z ,L X IL Z 

u 2u u 



2u 2 Au 2 z u 2 



+3a(B t B' t - fB z B' z + l -B t h\ + B z h\ - B t h x x - l -B t h z z ) J 



u=0 

(3.5) 
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Relevant terms of the Gibbons- Hawking term (12.251) are explicitly given by 



Scn ~ 32^ J j2^\u A ~ u KK + u KK 



u u u 

+—hlh x J + f -h\h z J + —h x x h z J 
u u u 



The counter term (12.261) also can be evaluated as 



^2 



+ ^ + ^ + ^[' (3 ' ?) 



up to 0(u , k ,uk) 



4 Pole Structures and Transport Coefficients from 
Hydrodynamics 



We now look for solutions of our set of equations fj3.4an - fl3.4dj) . and (I3.3al) and (13.3bl) . 



We will consider these equations of motion in low frequency limit so-called hydrodynamic 
regime. In this regime we could obtain the sound velocity, the diffusion constant for U(l) 
charge and the electrical conductivity from retarded Green functions. 

4.1 Master variables 

By using master variables derived by Kodama and Ishibashi in [19], the following field 
is introduced: 

* s 4 M 3/4 (4 ^-3/o (( 4 ''^ 2 - 3f ') K + 2/ K + "•'))■ <"> 
For the gauge field, the corresponding variable is given by 

A = 2a(^-h\ + 3hl-2B t '\ (4.2) 
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I6u 3 / 2 f(4b 2 k 2 - 3f ) 2 



In terms of these new variables $(w) and ^4(m), Einstein equations (I3.4ap - (13.4dl) and 
Maxwell equations (I3.3al) and (13.3b|) can be combined as follows: 

o = (u 1 ' 2 !®)' 

( -lQb 2 uj 2 u(4b 2 k 2 -3f) 2 

+f 2 (lQ(-b 4 k 4 + W8ab 2 k 2 u 2 + 162aV) 

+27f'{8b 2 k 2 + lQau 2 + 5/')) 
+Auf(4b 2 k 2 - 3f) (lQb 2 k 2 (b 2 k 2 + 3au 2 ) 

+f (8b 2 k 2 + 36au 2 + 9f )) |$ 

^2_ 3/02 { 3 /( 4&2fc2 + V + W ) +^(l6& 4 fc 4 - 9(f) 2 ) }A 

(4.3a) 

= (ufAj 



8mV4(452/ 



"——^(Atfk 2 - 3f) (u 2 - k 2 f) - ISauf^A 



f(4b 2 k 2 

-48au 3/4 f& 
4a 

+ u 1 /\Ab 2 k 2 -3f 



-^32b 4 k 4 u + Ylf(b 2 k 2 + 9aw 2 ) + 3f ( - 86 2 £; 2 w + 9/) J$. 

(4.3b) 



Next we would like to try to obtain decoupled equations from the equations (14.3aj) and 
( 14.3bj) . This will be done by introducing the following linear combinations of the variables: 

$ ± = a ± $ + PA, (4.4) 

where the coefficients a± and f3 are 

a± = C± — 3aw, 
u 1/4 
H 8 ' 

with the constants 

C± = (1 + a) ± v/(l + a) 2 + 4a6 2 fc 2 . 

As a result, we can obtain second order ordinary differential equations in term of these 
new variables, 

(u l/2 fY 

o = $± + ^7^r$± + v±*±, (4-5) 

U ' J 
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where potentials V±(u) are given by 

v - = Wf ^-3n { mWuiAb2k2 ~ 3/r 



-4uf{Ab 2 k 2 - 3f) (lQb 2 k 2 {b 2 k 2 - C±u + 3au 2 ) 

-Af(2b 2 k 2 + 3C±u - 9au 2 ) - 3(f) 2 
+f 2 {w(b 4 k 4 + 12C ± b 2 k 2 u 

-I08ab 2 k 2 u 2 + 54C±au 3 - 162aV) 

-24f'(b 2 k 2 -6C±u- 18aw 2 ) +9(f) 2 }}- 



(4.6) 



Considering the perturbative expansion with respect to small uj and k, it might be 
convenient to introduce new variables 

$ ± = H ± * ± , (4.7) 

where the factors F±(u) are 



H 



± 



m" 3 / 4 (for $ H 



— (for 



(1 + a) — -au 

so that the second order differential equations (14.51) are reduced to be much simpler forms 
to solve, 

= ^+ ( -^/f^ + V ± ^ ± , (4.8) 
where the potential V±(u,u, k) is newly defined from the original potential V±(u,u, k), 

V±(u, cu, k) = V±{u, cu, k) - V±(u, 0, 0). (4.9) 

4.2 Perturbative solutions 

Let us proceed to solve the differential equations. First we consider the equation for 
$ + (m). Following the usual way to solve differential equations, we impose a solution as 
= (1 — u) u F + (u) where F + (u) is a regular function at the horizon u = 1. Plugging 
this form into the equation of motion, one can fix the parameter v as v = ±iu/(47rT) 
where T is the temperature defined by eq. (|2.18|) . We here choose 

" = - ! Sr (4 ' 10) 
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as the incoming wave condition. We are now in the position to solve the equation of motion 
in the hydrodynamic regime. We start by introducing the following series expansion with 
respect to small uo and k: 

F+(u) = F +0 (u) + uF +1 (u) + k 2 G +1 (u) + u 2 F +2 (u) + 0(u\ ook 2 ), (4.11) 

where F +0 (u), F + i(u) and G + i(u) are determined by imposing suitable boundary condi- 
tions. The solution can be obtained recursiveljj*]. The result is as follows: 

F +0 {u) = C, (const.) (4.12a) 
^ / \ iCb f / 9 x 6Ki(u) 1 , 

= {** (i + « - - 2 ) - , (4.i2b) 

2 „ f 1 



G +1 («) = -CV -=^= - — — , (4.12c) 
3 tvl + 4a (l + a)w 

/" n Gh 2 
F +2 (u) = / d« 



(1 — u)(l + u — au 2 ) 

(1 - u)(l + au) log (l + u-au 2 ) 3(1 - u)(l + au)K 2 (0) 
" ~ U + 2(2 -a) 2 2(2 - a) 2 v / TT4^ 

1 + a)K 2 (l)w ^3 + (5 + 3a - 6a 2 + 2a?)u - 3au 2 ^jK 2 (u 



+ 



VI + 4a 2(2-a) 2 \A + 4a 



(4.12d) 



where 



KAu) = - logfl + u — au ) — log 1 



2 &v V l + v / T+4^ 

_ _ , , . / 1 + VI + 4a - 2aw 

K 2 (u) = log . 

V-l + VI + 4a + 2cra, 

Next, we shall study the equation of motion for Assuming again $_(tt) = 

(1— u) u F-(u) where F_(u) is a regular function at u = 1, the singularity might be extracted. 
We fix the constant as v = — iujji^nT) to use the incoming wave condition. We now impose 
a perturbative solution as 

F_{u) = F_o(u) + + k 2 G^(u) + u 2 F_ 2 {u) + O(oo 3 , uk 2 ), (4.13) 

and then we obtain the following resultf}]: 

F_ (u) = C, (const.) (4.14a) 



* The derivation of the solution is given in Appendix B. 
t The detail is given in Appendix C. 
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F-i(u) = {8(1 + a) 2 log(u) - (2 + a)(l + 4a) log (l+u- au 2 ) 



G_i(u) = C6 : 



2(2 

-2V1 + 4a(2 + 5a)ATi(u)}, (4.14b) 

2(l + a)(2 + a) log(u) 



2(l + a) 2 (l + a-^a«) (2 a) ' 

(l + a)(2 + a)log(l+ M -a M 2 ) ^ 2(2 + 5a + 6a 2 ) if ^ (4 ^ 



(2 -a) 2 

dM 2(2 - a) 4 (l + 4a) 3 / 2 (l - u)u (1 + u - au 2 ) 

x ^8(2 - a)(l + a) 2 (l + 4a) 3/2 w(l +u- au 2 ) log(w) 

-(2 - a)(l + 4a) 3/2 (4(1 + a) 2 + (2 - 3a - 8a 2 )w + (2 + 9a + 13a 2 )w 2 

-a(2 + a)(l + 4a)w 3 ) log (l + u - au 2 ^ 
+ (1 + 4a) 2 (2 + 5a)X 2 (0)(l - u) (4(1 + a) 2 + (2 - 3a - 8a 2 )w 

-(2 - a)aw 2 ) 

-2(1 + a) (2 - 2a + 41a 2 )fT 2 (l) ^2(1 + a) - 3aw) 2 

-(2 - a) (a(l + 4a) 2 (2 + 5a)w 3 - (2 - a)(l + 11a + 46a 2 + 18a 3 )w 2 
- (2 + 9a + 180a 2 + 224a 3 + 24a 4 ) u 

-4(1 + a) 2 (l - 10a - 2a 2 ))ff 2 ( M ) J. (4.14d) 

Using these, we can get behaviors for the solutions of $+(«) and around the 

boundary u — 0, 

H 1 -*^* 1 *^^*--)' (4 - 15a) 



p; iM f 1 fAi(l + a)bu 2(2 + a)b 2 k 2 

*- = c " 1/4 {tt^ + Hr^r - + J log(l,) + 

(4.15b) 

where the constant £)_ is 

D - ^ (2-a) 4 (l + 4a) 3 / 2 { " 27(2 " 

+4(l+4a) 3/2 (l + a) 3 log(2-a) 
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+4(2 - 2a + 41a 2 )(l + afK^l) 

Let us now consider the integration constants C and C . These could be estimated in 
terms of boundary values of the fields 

lim/i*(u) = (Mf, lim B t (u) = (B t )°, etc. 

Using equations of motion, the integration constants C and C are determined as 

C = 2{k2 1 _ 3uj2) {^ 2 (B t r + 3aku(B z r 

+ (1 + a) ( - k 2 {h\f + 2ku(h z t )° + (k 2 - u 2 ) {h x x f + u 2 (h z z )° 



(4.16a) 

(2-a) 2 abk(k(B t f + u(B z )°) 

C = 9 n ( m ( 4 - 16b ) 

2D p (u, k) 

where 

D p (u, k) = 2(2 + a)bk 2 - 4z(l + a)u - (2 - afbD^u 2 . 

In the equations (14. 16a,}) and ( I4.16bl) . one can observe the existence of the sound and 
diffusion poles in the complex tu-plane. 

4.3 Retarded Green functions 

Let us evaluate the Minkowskian correlators. The relevant action is given by the sum of 
three parts ([33]), flM]) and fl3TTj) . 

S = So + SqU + Set 

- _J [ d k S-h z h z/ + -Z? 3 ^' 4- -h l h x ' + —h l h z ' 

+ l h:h f + L KK > + L h z x' + i-h z x 

u u 2u u 

z\2 



Au 2 
3 



(/ - v 7 /) (*4) 2 + 7^ (3/ - «/' - 3 a/7) 
(/ - V 7 /) (^) 2 - ^ (6/ - uf - 6y/f) 



u 2 f 

(e/ - «/' - 6x77) ^ - ^ (3/ - uf - 3 v^7) kk 



+3a(B t B' t - fB z B' z + X -B t h\ + £^ " ^ " } 



u=0 

(4.17) 
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Using equations of motion and solutions of and derivatives of h's and -B's 

can be expressed in terms of their boundary values: 

^^ra{- M (* (B ' )0+ - (B "»°) 

+(l + o.) (V(A;) - 2uk(Kf - w 2 (2(ft*)° + (K)°)) }, (4.18a) 



+(1 + a)[- fc 2 (ft|)° + 2kMhir + w a (2(fc3° + («)") )l (1.1 81,) 
+(1 + a) ( - fc 2 (ftD° + 2to W)° + " 2 (2(/»3° + (Kf)) } , (4.18c) 

^'-tct{- 9 ^(^ )0+ ^ )0 ) 



-(1 + a)k(2k(hff + W (-3(fci)° + 2(Kf + ('»;)")) }, (4.18d) 

B '- 2(l+a)(L3^) {- 9a K^ + ^>°) 

+ (1 + a) ((2A; 2 + 3cu 2 )(^)° - 6k^t)° 

-3^ 2 (2(^)° + (^) )) 

(2-a) 2 &^(£ t ) + ^^ ) 

(A 18e) 

(l + a)D p (u,k) ' { } 

^- 2(l + fl )(P-3^ { 9 -(^ + ^°) 



+(1 + a) ( - 3ku(^)° + 2& 2 (/if) + kuj{2{h%f + (^)°~ 



(2-a) 2 cWA;(5 t ) + u;(5 2 ) 



+ (1 + a)D p (uj, k) ' (4 ' 18f) 

Substituting these expressions to the surface term (14.1 7ft . we can read off the Green func- 
tions defined by (12.51) . Through the counter terms, the singularities around the boundary 
vanish completely. The results are listed below in Table [U [2] and [3l 

In the final expression we rescaled the gauge field (-B„)° to the original one (A^) = 



;2 



(Bfj) and raised and lowered the indices by using the flat Minkowski metric rj^ 



Z 1 
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Table 1: G,» ^ / ( ^gggg^ 3^) ) ' 





tt 


XX 


zz 


tz 


tt 


3 (5A; 2 - 3lo 2 ) 


6(k 2 + co 2 ) 


3 (k 2 + uo 2 ) 


24(k 2 + to 2 ) 


XX 




IQlu 2 


2 (fc 2 + to 2 ) 


16ku 


zz 






-k 2 + 7lu 2 


8kuj 


tz 








4(k 2 + 9a; 2 ) 



Table 2: G** * / 3 ^ 2) 







XX 


zz 


tz 




3A; 2 


2k 2 


k 2 


6kcu 




3ku 


2ku 


ku 


6uj 2 



Table 3: °* * 1 (4e 2 (l + a)6 2 (fc 2 - fr^ + ^ p( J^) ) ) ■ 





t 


z 




k 2 


ku 


z 




cu 2 
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diag(— , +, +, +) in four-dimensional boundary theory. Taking the limit which the charge 
goes to zero, the correlators for energy-momentum tensors coincide with the known ones 
in [15]. In this limit, the correlators for the energy-momentum tensor and the U(l) current 
vanish, while ones for the U(l) currents have no sound poles, as we could see in the case 
of vector type perturbation for (1,0,0) R-charged [22] and RN-AdS black holes [20]. 

Remember 1 3 /k 2 = N 2 /(Att 2 ) and it should be noticed that the factor l/e 2 = N 2 /(16tt 2 ) 
for the R-charge since e = 2k/ I in that case, while l/e 2 = N c Nf/(4ir 2 ) for the brane 
charge [29]. 

4.4 Transport coefficients 

From the obtained Green functions, we can observe the value of the speed of sound without 
the medium effect, 

v. = -jg. (4,9) 

One should notice that there is no effect of the charge on the sound velocity. 

We can also find the diffusion pole in the current- current correlators. The diffusion 
constant can be read off 

Da = (4-20) 
A 2(1 + a) K J 

It should be compared with the diffusion constant for gravitation fields D H obtained in the 
previous work [20], 

so that the relation between them is 

D A = (2 + a)D H . (4.21) 

In the chargeless case, we can reproduce the result in [15]. 

The electrical conductivity a of the medium could be also determined by the current- 
current correlators via Kubo formula, 

2 

a = - lim p-Tmf&'GiAu, k — 0)\ 

where e-^ is a four-dimensional gauge coupling. Together with the result for vector type 
perturbation [20], we can obtain 

Ze|(2 - a) 2 



24e 2 (l + a) 2 6 
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x 3 



e 2 J 2(1 + a) 2 

We can also access to the charge susceptibility 5 defined by 



l\^(2-a) T 



(4.23) 



T (volume) 

Using Green function Gt t{k) which might give an expectation value of Q 2 , one can obtain 
the following relation in thermal equilibrium, 

(^i)=/l(- Im ( G '^^°»))"*»' < 4 - 24 ' 

where n^oj) = e u,/r_ l is Bose-Einstein distribution function [31]. From Table 3, we can 

see 

r2;^2 / ,,n.i,2 



Tm(G t t (uj, k) 



A-nHT 2 ( uoD A k 



\ v ' V e 2 (l + a)(2 + a)\uj 2 + (D A k 2 ) 2 , 

It should be noted that a quantity D A k 2 / (uj 2 + (Z^A; 2 ) 2 ^ approaches to 2n5{u) for — > 
limit. Taking the relation (14.231) into account, we can read off the charge susceptibility H 

as 

Z \ 4vr 2 2 . 

^7 7 w -T . 4.25 

3 2 ; (l + a)(2 + a) 1 ; 

We can then observe that Einstein relation 

a/(e E 2 E) = D A (4.26) 

holds exactly. (See also [32] .j^ 

In R-charge case, taking the charge-free limit, the electric conductivity and the charge 
susceptibility coincide with the results in [31]. 

It is interesting to express physical constants in terms of the boundary variables: the 
temperature and the chemical potential. In fact, it is easy to verify that 

« = 2 ; 4 &=(-L] 1 (4.27) 

1 + v/l+4(/VT) 2 ' \kTJ 1 + v /l+4(/ t /T) 2 ' 

K 1 

where ft = —= — fi. Notice that for the R-charge, ft = — -j=- //, while for the brane 
y 37re/ 2v37r 

/ N f 

charge, ft = \ -—- — ^u. The behaviors of the diffusion constants D A and Dh, the electrical 

V oNrlt 



* It is interesting to notice that (dp/d/j,)T gives a different value for S given in eq. (|4.25[) . The authors 
thank J. Mas and J. Shock for pointing this out. 
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conductivity a and the charge susceptibility S are drawn as functions of the ji/T in Figured], 
El [3] and H] respectively. Notice that for the fixed temperature, all of them are decreasing 
functions of the chemical potential. One should notice that there is no upper bound of 
ji/T for any of these quantities unlike (1,0,0) charged black hole studied in [22]. 




Figure 3: a/{Tle\/e 2 ) vs. ji/T. Fi S ure 4: S /( T ^/e 2 ) vs. ji/T : Notice the 

rapid change between two finite values as 
T runs from to oo 

It is particularly interesting to observe that the charge susceptibility modulo T 2 factor, 
which is an indicator of the degree of freedom, shows rapid change between low and high 
temperature (density) indicating a mild phase transition. Such behavior does not exist for 
chargeless case. See also [33]. 
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5 Conclusions and Discussions 



In this paper, we worked out the decoupling of scalar modes of the charged AdS black 
hole background in 5*0(2) basis for the mode classifications. We also perform the hydro- 
dynamic analysis for the holographic Quark-Gluon Plasma system. Master equations for 
the decoupled modes are worked out explicitly. The sound velocity is not modified by the 
presence of the charge. We calculated the diffusion constants, the charge susceptibility and 
the conductivity as a consequence and observed that Einstein relation holds between them. 
These transport coefficients are modified due to the charge effect. Interestingly, the sus- 
ceptibility modulo T 2 factor, which is an indicator of the degree of freedom, shows rapid 
change between low and high temperature (density) indicating a mild phase transition. 
Such behavior does not exist for chargeless case. 

One can give an explanation of hydrodynamic mode in meson physics. In our interpre- 
tation, the Maxwell fields are the fluctuations of bulk-filling branes, therefore they should 
be interpreted as master fields of the mesons. Then hydrodynamic modes are lowest lying 
massless meson spectrum. In terms of brane embedding picture, this massless-ness is due 
to the touching of the brane on the black hole horizon. Near the horizon, the tension of 
the brane is zero due to the metric factor and it can lead to the massless fluctuation. Then 
the massless spectrum can not go far from the horizon in radial direction. In this picture, 
hydrodynamic nature is closely related to the near horizon behavior of the branes. We 
will discuss the spectrum of meson mode by considering the quasinormal mode [34] of the 
vector modes in elsewhere. 
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Appendix A Results for the vector and the tensor 

type perturbations 

Appendix A.l Vector type perturbation 

In the vector type perturbation, independent variables are 

h t Jx), h zx (x), AJx). 



One can observe the diffusion constant for the metric perturbation Dh as 

b 



D 



H 



(A.l) 



2(1 + a) 

We list the retarded Green functions in Table HJ [5]and[6j Notice that 1 3 /k 2 = N 2 /(4ir 2 ) 



Table 4: G„ / 



l6K 2 b 3 (ico - D H k 2 ) 





tx 


zx 


tx 


k 2 


—ujk 


zx 




CO 2 



Table 5: G« * / 



Z//, 



4e 2 b 2 (ico - D H k 2 ) 





tx 


zx 


X 


—Hlo 


b I 

-uuk 

(1 + a) 



/ I f 3a (2-a) 2 b\\ 

Table 6: G* * / ^ 4e2(1 + o)&2 [ iuJ _ p Hk 2 ~ 2(1 + fl) ))■ 





X 


X 


iuj 



and l/e 2 = N 2 /(16n 2 ) for the R-charge, l/e 2 = N c Nf/(4:7r 2 ) for the brane charge. From 
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the Green function for U(l) currents G x x (u, k), one can read off the thermal conductivity 
kt via Kubo formula, 

«r = lim -Wc. jfc = 0)) = 2vr 2 ^^. (A.2) 

Appendix A.2 Tensor type perturbation 

In the tensor type perturbation, an independent variable is just 

By using Kubo formula, one can obtain the shear viscosity rj as 

1 / ^ I 3 

= — lim — Im( G r 



V 



•xy xy 



where the retarded Green function is given by 

Table 7: G** / 



(u,k = 0) 



16k 2 6 3 ' 









icj + bk 2 



(A.3) 



One can confirm the universal (within Einstein gravity!*] ) ratio that is the ratio of the shear 
viscosity to the entropy density s, 

(A A) 



J] 1 

s An 



Appendix B Perturbative Solutions for <l> + 

Substituting the equation (14.111) into the equation (14.81) . one can read off one for F +0 (u), 

= (V x (l -u){l + u- au 2 )F' + ^j . (B.r 
A general solution is given by 



F +0 (u) =C + £) 1 2 log (l-uj - log (l + u- au< 



3K 2 (u) 
VI +4a 



:b.2) 



* If one consider higher derivative corrections to Einstein gravity, the viscosity bound could be modi- 
fied [35-39]. 
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Constants of integration Cq and D should be determined to be a regular function at the 
horizon. So we here choose Dq = and set 



F +0 (u) = C = C. (const.) (B.3) 
By using this solution, one can get an equation for F + i(u) from the equation (14. 8p . 

= (u-^l -«)(! + «- au 2 )F' +1 ) - i C ^_ + a ^ ] - (B.4) 
A general solution is 

i(c(l - a)b - (2 - a)D^j log (l - u) 



F +1 {u) = d + 



(2 -a) 



2 



i (Cb + (2 - a^Dj) log ( 1 + u - an 2 ] 3i (cb + (2 - a)Dx) K 2 (u) 

' ) + — ^ '— . (B.5) 



2(2 -a) 2 2(2 - a) 2 VTTia 

Removing the singularity at the horizon, the constant D\ should be 

D x = Cb- ~ ° 



2- a 

We also impose a boundary condition F + %(u = 0) = so as to fix the constant C\. 
Therefore the final form is 

^-^H^-^-MY (B6) 

A differential equation for G + \{u) is 

, , Cb 2 U- (l + a)u 2 ) 
0=(u-\l-u)(l + u~au 2 )G' +1 ) + V 3(1 + qK J . (B.7) 

A general solution is 

G +1 [u) = C x 



4C(2 - a)b 2 



6(2-a)(l + a) t u 

(C(l + 5a - 2a 2 )b 2 - 9(1 + a)5 x ) K 2 {u) 
+~ / " 

+ (C(l - a)b 2 + 3(1 + a)5i) 

x^21og(l-u) -log(l + u-au 2 ^ j (B.8) 
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and the constant D 1 can be fixed as 

C(l-a)b 2 



D 1 



3(1 + a) 



From the condition (uG + i)'\ u=0 = 0, we can fix the constant C\. The final form of the 
solution becomes 

2 , f K x {u) 1 



G +1 (u) = -Cb 2 { - - } ■ (B.9) 

+U ; 3 Iv^T^ {l + a)u ( 



A differential equation for F +2 {u) is 



= 



u 



- u) (1 + u - au 2 ) 



Cb 2 log (l - uj ibF+iju) ib\og{\ - u)F' +1 {u) 
(2-a) 2 (l-w) (2-a)(l-u) 2-a 
Cb 2 



u 2 {l — u)(l + u — au 2 ) 
Integrating over u, we have 



(B.10) 



K2H 



2(2 - a) VI + 4a(l - u)(l + 

x |2 V / TT4^(2 - a) 2 (Cb 2 - D 2 u) 
-3Cb 2 K 2 (0) (l + u-au 2 ) 

+Cb 2 (3 + (5 + 3a - 6a 2 + 2a 3 )u - 3aw 2 ) ir 2 (u) 

+ 4a6 2 (l - u)(l + aw) log (l + « - aw 2 ) J (B.ll) 

and the constant -D 2 can be fixed as 

I 2(2-a) v / TT4^ V / TT4^ J 
so that ((1 - u)F| 2 (u)) | u=1 = 0. Then, we find 

F| 2 (u) = Co 2 + £>(«). (B.12) 

Hence, we can write the solution of F +2 (u) as 

Cb 2 



'+2 



d-u 



(1 — u)(l + u — au 2 ) 
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(1 —u)(l + au) log ( 1 + u — au 2 ) 

x<! l-u+ ; ^ - 

2(2 -a) 2 

3K 2 (0)(l-u)(l + au) (l + a)K 2 (l)u 



2(2 - a) 2 Vl + 4a VI + 4a 

(^3 + (5 + 3a - 6a 2 + 2a 3 )w - 3au 2 ^)K 2 (u) 
+ 2(2 - a) 2 v / T^1a 

which satisfies a boundary condition -F +2 (0) = 0. 

Appendix C Perturbative Solutions for <£_ 

For F_ (u), one can get an equation 

= (u(l -u)(l+u- au 2 ) (l + a- ^auj F^^j ■ 
A general solution is given by 

F- {u) = Co 

+ A) I (2 - a) log (l - u) - 4(1 + a) 2 log (u) 

l, n , Nl /. 2 \ VI +4a(2 + 5a)i\r 2 (M) 
+-(2 + a)(l + 4a) log [1 + u - au 2 ) - - y — ; v ; 

Since the function F_o(w) should be regular at the horizon, we choose Dq = and 

F_ (u) = C = C. (const.) 
Substituting the solution to the equation (14.81) . we get an equation for F_i(u), 

= (u(l -u)(l + u- au 2 ) (l + a - \u) * F'_ \ 

iCb(2(l + a) + (4 + 7a)u - 6a(l + a)u 2 + 3aV 



/ 3 \ 3 

2(2 - a)[l + a- -au) 

A general solution is given as 
F-i(w) = Ci 

27 



54(2 - a)a 2 

x |8(1 + a) 2 (d(l + Aa)b - 27L>i(2 - a)a 2 ) log(u) 

-2 (c(2 + 7a + 23a 2 )6 - 27D 1 (2 - a) 2 a 2 ) log (l - u) 

-(2 + a)(l + 4a) (c(l + 4a)6 - 27L>i (2 - a)a 2 ) log (l + u - aw 2 ) 

+v / TT4a"(2 + 5a) (d(l + 4a)6 - 27D 1 (2 - a)a 2 ) if 2 («)|. (C.5) 

The constant of integration D x should be 

_ Cb{2 + 7a + 23a 2 ) 
1 ~ 27(2 - a) 2 a 2 ' 

so that the singularity at the horizon could be removed. In addition, we require the 
condition 

F.xH-logH lim 

^ log(u) 

to fix the constant C\. Then we get the final form of the solution 



= 

u=0 



F-M = 2( 2_ fl)2 {8(l + a) 2 log(u) - (2 + a)(l + 4a) log (l + u - aw 2 ) 

-2 v / TT4^(2 + 5a)iri(M)|. (C.6) 

Similarly we have a differential equation for G-i(u), 

= (^u{\ - u){l + u - au 2 )(l + a - ^auj 2 G'_^j 

Cb 2 { 
+ = — -g <^ - 4(1 + a) (2 + 6a + 3a 2 + 2a 3 ) 

8(1 + a) (l + a - -auj ^ 

+2a(10 + 30a + 57a 2 + 10a 3 )w 

-42a 2 (l + a)V + 21a 3 (l + a)« 3 j. (C.7) 

A general solution of this equation can be obtained 



G-i(u) = Cl Cah 3 

a(l + a) \ \ + a — -^ a j 

1 



54a(l + a)(2-a) 
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x |4(1 + a) 2 (2 - a) {iCb 2 + 54D 1 a(l + a)) log(u) 

- (2C(U + 13a + 17a 2 )o 2 + 54Z)i(2 - a) 2 a(l + a)) log (l - uj 
-(2 + a) (d(7 + 11a - 14a 2 )6 2 + 27Z>ia(2 - a)(l + a)(l + 4a) 
x log ^1 + u - au 2 j 
1= f (7(14 + 57a + 81a 2 - 16a 3 )6 2 



VI + 4a V 



+27Dia(2 - a)(l + a)(l + 4a) (2 + 5af)K 2 (u) 1. (C.8) 



The constant of integration D\ might be fixed to remove the singularity at u — 1, 

~ Cb 2 {U + 13a + 17a 2 ) 
1 ~ 27a(l + a)(2-a) 2 

Another constant of integration C\ is fixed to satisfy the condition 



G-i(u) — log(-u) lim 
The final result of the solution is 



u=0 



G-i(u) = -Cb 2 



9a 2 u 



o y log(-u) 

6(1 + a)(2 + a) 



0. 



2(1 + a) 2 (l + a- -au) 
3(l + a)(2 + a) 



(2 - aY 



log(u) 



(2 -ay 

A differential equation for F_ 2 (u) is 
u(l + u — au 2 ) 



log (l + u — au 2 ) + 



6(2 + 5a + 6a 2 
{2-a) 2 y/T+ta 



Mu) . (C.9) 







/ 3 \ 2 

(2 - a) 2 (l + a - -aw J 

x ^(2 - a) 2 (l - u)F'_ 2 {u) + i(2 - a)6F_!(u) 

-i(2 - a)6(l - u) log (l - u^Fi^u) + Cb 2 log (l - u) 



Cb 2 



^1 + a — 2 aM ) (1 — m) (l + -a — aw 2 ) 



(CIO) 



Integrating over u, we have 

3 N 2 



1 + a 



6 Y 

--au) 



(1 — u)u (1 + u — au 2 ) 
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f „ 18Cab 2 

X \ D2 + 7 3~T 

L {2-a) 2 (l + Aa)(l + a--au} 



4C(1 + a) 2 b 2 u(l + u- au 2 ) log(u) 

3 2 

(2 - af(l + a - -au) 

/ ... . . . n \ 



2Cb 2 



'{2-af 



(2 + a)(l + Aa)u(l + u - au 

_V *{l + a--au) j 

CVma(2 + 5a)b 2 K 2 (0)u(l + u- au 2 ) 

2( y 2-a) 3 ^l + a-^a U y 

2Cb 2 ( 
+ (2-a)3(l + 4a)3/ 2 ( 1 - 2a(5 + a) 



log ^1 + u — au 2 ^ 



(2 + 5a)(l + 4a) 2 w(l + u - au 2 )\ T \ 

/ 3 — \2 J K ^ U ) f> 

4\l + a--au) / > 

(C.11) 

and the constant D 2 can be fixed as 

° 2 = (2-ay(l + Aa)y 2 \ VTT ^V {1 + 4 ° )(1 + ^ l0g(2 " a) " 9 ° (2 " a) ) 

+ ^(2 + 5o)(l + 4a) 2 K 2 (0) - (1 + a) (2 - 2a + 41a 2 )X 2 (l)|, 

so that ((1 — u)F'_ 2 (uj) | = 0. Then, we can write the solution as 

/U Qff 
dM 2(2 - a) 4 (l + 4a) 3 / 2 (l -u)u(l + u- au 2 ) 

x ^8(2 - a)(l + a) 2 (l + 4a) 3/2 w(l + u - aw 2 ) log(u) 

+ (2 - a)(l + 4a) 3/2 (a(2 + a)(l + 4a)w 3 - (2 + 9a + 13a 2 )w 2 

-(2 - 3a - 8a 2 )w - 4(1 + a) 2 ) log (l + u - aw 2 ) 
+ (1 + 4a) 2 (2 + 5a)K 2 (0)(l - u) + a) 2 + (2 - 3a - 8a 2 )w 

— a(2 — a)u' 
-8(1 + a)(2 - 2a + 41a 2 )iT 2 (l) (l + a - 

-(2 - a) (a(l + 4a) 2 (2 + 5a)w 3 - (2 - a)(l + 11a + 46a 2 + 18a 3 )w 2 



, 2 

-au 
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-a(l + 11a + 46a 2 + 18a s )u 

-4(1 + a) 2 (l - 10a - 2a 2 ))K 2 (w) j. (C.12) 
Expanding this expression around u = 0, we have 

F_ 2 (u) = Cb 2 D4l + a) log(u) + C 2 + (C.13) 
where C 2 is an integration constant. 
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